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[1] Lineare Gleichungssysteme losen mit DERIVE for Windows, B Kutzler, bk-teachware 

[2] Lineare Gleichungssysteme losen mit dem TI-92, B Kutzler, bk-teachware 

[3] Optimierungsaufgaben, graphisch, numerisch und analytisch mit dem TI- 92, J Bohm, bk-teachware 

[4] Solving Systems of Linear Equations with DERIVE for Windows, B Kutzler, bk-teachware 

[5] Solving Systems of Linear Equations with the TI- 92, B Kutzler, bk-teachware 


bk-teachware, Softwarepark, A-4232 Hagenberg, Austria, email: bk.teachware@swp.co.at 
Chartwell-Yorke, 114 High Street, Belmont, Bolton, UK, email: Philip.Yorke@chartwellyorke.com 


AUTHORS SOUGHT 

Authors are sought for the bk teachware series: " Support in Learning”. Each text is written for the average, 
technology inexperienced teacher, explaining how to use one of the following tools for teaching: TI-92, TI-89, 
CBL, CBR, DERIVE, CABRI GEOMETRE II. 

We offer: 

- quick publication, professional production, attractive sales prices 

- 15% of net income as royalty 

- publication with ISBN number 

- experienced editors: Brian Denton (Editor-in-Chief) Terence Etchells 
(Assistant Editor) Philip Yorke (Desk Editor) 

- marketing through distributors in the UK, USA, Scandinavia, and Benelux 

- option for translating and marketing in other language regions (German, French, Italian, Japanese, ...) 

The booklets of the bk teachware series have a format of 17 x 24 cm and comprise 32-64 pages. A diskette can 
be attached to the text. They are priced between 4 and 10 Pounds. 

One booklet should treat one topic, ideally should be self-contained, and should describe a lesson in all details. 
The following booklets can serve as models for your own:- 

SL-03: Solving Linear Equations with the TI-92 

SL-04: Solving Systems of Linear Equations with the TI-92 

SL-05: Solving Linear Equations with Derive for Windows 

SL-06: Solving Systems of Linear Equations with Derive for Windows 

We would be pleased to send you sample pages from one of these booklets. Please contact the desk editor: 

Philip Yorke: Philip.Yorke@ChartwellYorke.com 
Please send manuscripts to the Editor-in-Chief: - 
Brian Denton: B.H.Denton@livjm.ac.uk 

Yours sincerely 

Bernhard Kutzler 
Publisher 




Interesting Web Sites: (only the red ones are still valid!) 

http://cs.gettyburg.edu/~leinbach/DRVTI_Conf/ConfPrg.html Three keynotes from the Gettysburg 

Conference 

http://hutchinson.belmont.ma.us/tth/tth.html A tool to translate TEX into HTML (free!!) 
http://www.imaxx.net/~gdorner/ G Dorner's pdf-files concerning a thrilling TI-92 book 

http://www.brunel.ac.uk/~mastmmg/ssguide/sshome.htm An Online Study Skills guide 
http://www.bham.ac.uk/ctimath Always interesting 

http://www.msri.org/realvideo/ 100 lectures (recommended by cti Birmingham) 

http://www.ti.com/calc/docs/t3ww The T 3 site (Teachers Teaching with Technology) 
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Liebe DUG Mitglieder, 


Dear DERIVE and Tl-friends, 


In der Nummer 1 nach Gettysburg mochte ich 
mich vorerst auch an dieser Stelle im Namen alter 
Teilnehmer bei den Organisatoren der diesjahrigen 
Konferenz fur ihren Einsatz recht herz-lich bedanken. 
Pat und Carl Leinbach haben es verstanden, die Tage 
vom 15.- IJ.Juli zu einem Fest fur die DERIVE & TI- 
92 Familie zu machen. Beth und ChristyAnn haben 
versucht, jeden Wunsch zu erfullen (sogar den absur- 
den Wunsch nach 7 Tennisschlagernl! 1^7. "You 
crazy Europeans with your 1 ", Carl L.). Leider konnte 
Bert Waits nicht bis zum Ende bleiben - er versaumte 
ein groflartiges Picknick auf Old Mac Leinbachs 
Farm. Wir haben ihn einige Zeit spater in seinem "Pa- 
rad ies" auf Seabrook Island besuchen durfen - noch 
vor dem Wirbel sturm Bonnie. (Bert, dein Sand dollar 
ist noch immer ganz.) 

In dieser Nummer ist noch ein Beitrag vom 1997 
Symposium in Sard zu finden: Terence Etchells zieht 
auf spitzfindige Weise aus numerischen Daten 
Schlusse auf deren logischen Zusammenhang. Ich 
habe ganz bewusst den Abschluss von D. Halprins 
Glove Osculants fur die nachste Ausgabe aufgeho- 
i ben, denn ich mochte, dass Sie seinen Beitrag be- 
quem anhand der files auf der Diskette des Jahres 
1998 nachvollziehen konnen. Aus dem gleichen 
Grund habe ich einen geplanten Artikel von S. Abra- 
ham Ibrahim uber "Modified Bessel Functions" in den 
DNL#32 verschoben. Es wird auch noch einen Bei- 
trag von S. Welke geben, der zeigt, wie sich DERIVE- 
Objekte auch ohne ACROSPIN animieren lassen. 

Ich mochte Sie hier auf die Moglichkeit hinwei- 
sen, eine preiswerte CD-ROM mit den Vortragen und 
Workshops der 3.lnternationalen DERIVE und TI-92 
Konferenz bei Jerry Glynns MathWare zu beziehen 
(Seite 8). 

Viele liebe Grufle 

Josef Bohm 



This is issue #1 after Gettysburg, and on behalf 
of all delegates I want to express my thanks to the 
organizers for this tremendous event. Thanks to Pat 
and Carl Leinbach, the days from 15-17 July be- 
came a festival for the DERIVE & TI 92 family. Beth 
and ChristyAnn actually fullfilled each and every 
wish (even the absurd request for 7 tennis rackets! ! 
1^7. ’’You crazy Europeans with your 1", Carl L.). 
Unfortunately Bert Waits could not stay until the end 
- he missed a great Picknick at Old Me LeinbaclTs 
Farm. We visited him and his wife in their "Paradise" 
on Seabrook Island - long before Bonnie, the hurri- 
cane. (Bert, your sanddollar is still unhurt.) 

In this issue you will find the final contribution 
from the 1997 Symposium in Sard. Terence Etchells 
draws logical conclusions from numerical data in a 
very inspiring and tricky way. I leave the last part of 
D. Halprin's Glove Osculants for the next issue, be- 
cause I think having the Diskette of the Year 1998 in 
your hands - or drive -, it will be much easier for you 
to reproduce his ideas. For the same reason I have 
shifted a planned article from S. Abraham Ibrahim 
about "Modified Bessel Functions" to DNL#32. 
There will also be a nice contribution from S. Welke, 
showing how to animate DERIVE objects even with- 
out ACROSPIN. 

Ed like to point out that there is opportunity to 
purchase an inexpensive CD-ROM containing the 
lectures and workshops of the 3rd International 
DERIVE and TI-92 Conference from Jerry Glynn's 
MathWare (page 8). 

With best regards 

Josef Bohm 



CSME f 99 

CALL FOR PAPERS 

Computer-Supported Mathematical Education 
(Didactic, Mathematical, and Software Technologic Aspects) 


August 23-25, 1999 
Castle of Hagenberg, Austria 


The submissions for CSME'99 are handled preferably via web or email. To use the web, carefully fill out the form available 
from the site of the conference: 

http://www.risc.uni-linz.ac.at/conferences/summer99/csme99/ 

Information: Editor of the DNL, CSME99-submit@risc.uni-linz.ac.at, hheugl@netway.at 
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The DERIVE-NEWSLETTER is the Bulle- 
tin of the DERIVE User Group. It is pub- 
lished at least four times a year with a con- 
tents of 44 pages minimum. The goals of 
the DNL are to enable the exchange of ex- 
periences made with DERIVE as well as to 
create a group to discuss the possibilities 
of new methodical and didactical manners 
in teaching mathematics. 

We have established a section dealing with 
the use of the 77-92 and we try to combine 
these modem technologies. 


Editor: Mag. Josef Bohm 
A-3042 Wiirmla 
D'Lust 1 
Austria 

Phone: 43-(0)660 3136 365 
email: nojo.boehm@pgv.at 


email: josef.boehm@blackboard.at 


Contributions: 

Please send all contributions to the Editor. 
Non-English speakers are encouraged to 
write their contributions in English to rein- 
force the international touch of the DNL. It 
must be said, though, that non-English ar- 
ticles will be warmly welcomed nonethe- 
less. Your contributions will be edited but 
not assessed. By submitting articles the au- 
thor gives his consent for reprinting it in 
DNL. The more contributions you will 
send, the more lively and richer in contents 
the DERIVE Newsletter will be. 


Next issue: December 1998 

Deadline: 15 November 1998 


Preview: Contributions for the next issues 

3D-Geometry, Reichel, AUT 

Linear Programming, Various Projections, Word Problems, Bohm, AUT 

A Utility file for complex dynamic systems, Lechner, AUT 

Examples for Statistics, Roeloffs, NL 

Fractals and other Graphics, Koth, AUT 

Implicit Multivalue Bivariate Function 3D Plots, Biryukov, RUS 

Riemann, a package for the TI-92, Bohm & Propper, AUT/GER 

Quaternion Algebra, Sirota, RUS 

Parametric 3D-Plots with DERIVE and 3DV, Welke, GER 

Various Training Programs for Students on the TI-92, Bohm, AUT 

A Critical Comment on the "Delayed Assignation" := =, Kurnmel, GER 

Sand Dunes, River Meander and Elastica, The lighter side Halprin, AUS 

Type checking, finite continued fractions, , Welke, GER 

Evaluation of the Modified Bessel Functions, Ibrahim, ESP 

On the Resolution of the Linear Differential Equation Candel, ESP 

TI-92 Statistics with DERIVE, Bohm, AUT 


and 

Setif, FRA; Vermeylen, BEL; Leinbach, USA; Weth, GER; Kirmse, GER 
Aue, GER; Koller, AUT; Mitic, UK; Schorn, GER; Santonja, ESP; 
Dorner, USA, Propper, GER and and and 


Impressum: 

Medieninhaber: DERIVE User Group, A-3042 Wiirmla, D'Lust 1, AUSTRIA 
Richtung: Fachzeitschrift 

Herausgeber: Mag.Josef Bohm, Herstellung: Selbstverlag 








D 


-N-L#31 


G P Speck: Concentric Curve Shading 


P 3 


CONCENTRIC CURVE SHADING 

G P Speck, Wanganui, New Zealand 

In the article ’’Broken Lines and Shaded Areas” in DNL#17 on pp 16-18, Josef Bohm gave a very 
simple elegant SHADE ( ) function which can be used to shade a region with line segments parallel to 
they-axis, where the region is bounded by curves of the formy = F(x) and y = G(x). Similarly, a triv- 
ial modification of the SHADE ( ) function gives a function which can be used to shade a region with 
line segments parallel to the x-axis, where the region is bounded by curves of the form x = F(y) and x 
= G(y). 

Of course, we can make the transparent observation that if a few marks (such as x’s) in a region are an 
adequate pointer to the region for some particular purpose, then these few marks can be placed in the 
region easily by devices such as the ( G r I ) + ( A) annotation device of DERIVE for Windows ( DfW ). 

In this article, shading regions bounded by curves given parametrically with concentric curves is con- 
sidered both as a way of creating a rather artistic pointer to a region, as well as a topic of interest in 
and of itself Of course, curves given parametrically are, in general, much more complex than curves 
given in the form y = F(x) or x = G(y). And even though it may be theoretically possible (in some 
cases) to reduce a given parametrization to a finite number of curves of the formy = F(x) or y = G(x), 
it may be EXTREMELY tedious to try to actually produce such a reduction, even with the aid of 
software such as DERIVE. 

Herein we will consider five fairly elementary parametrizations in SHADING l.MTH through 
SHADING5.MTH. These MTH file names will serve both as our names for our parametric curves, 
which will be given first, as well as the names for the files which show how the curves with shading 
were constructed, which we will give together with the graphs. 

The regions bounded by parametric curves in SHADING1.MTH, SHADING2.MTH and 
SHADING3.MTH are convex sets, for which our method of concentric curve shading will be seen to 
be appropriate since any ray from an interior point of a convex set intercepts the boundary in one and 
only one point. Thus, though the regions which we wish to shade in SHADING4.MTH and 
SHADING5.MTH are not convex, our method of concentric shading will still be seen to be appropri- 
ate since, in each case, there exists a particular point interior to the non-convex region such that each 
ray from the point intercepts the boundary in one and only one point. 

Also, since my original version of DfW will not PLOT curves given via I F ( ) statements (a real 

shortcoming!), I have defined characteristic functions which enable me to define piecewise parametri- 
zations without the use of IF ( ) statements. (Incidentally, I have found these characteristic functions 
to be very useful elsewhere - e.g., in shoring up DERIVE shortcomings in limit evaluations.) 


My self-explanatory (hopefully) examples follow. 
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#1: SHAD_UT . MTH 


Method via Characteristic Functions for my version of DfW 
which will not plot functions given via IFQ statements 


#3: CHARACTERISTIC FUNCTIONS 


#4: <x(x) := FLG0R(x + l)-x(0, x + 1, «0 - FLOOR(x) . x (0, x, «0 

#5: CHPOINTCx, a) := a(x - a) -a (a - x) 

#6: CHGRTEQ (x , a) := a(x - a) 

#7: CHGRT(x, a) := a(x - a) - CHP0INT(x , a) 

#8: CHLESS(x, a) := 1 - a(x - a) 

#9: CHUNEQ (x , a) :=CHLESS(x, a) +CHGRT(x, a) 

#10: CHLESSEQ (x , a) := 1 - a(x - a) + CHPOINTCx, a) 

#11: CH_QLCL(x) := a(x)-a(l - x) 


# 12 : 


#13: 

#14: 

#15: 


CHCL0SED(x, a 

CH0PENRT(x, a 
CH0PENLF(x, a 
CH0PEN(x, a, 


, Id) := CH_01_CL 

k b - 

, b) := CHCL0SED(x, a 
, b) := CHCL0SED(x, a 
b) := CHCL0SED(x, a, 


(x - a) 

a j 

, b) - CHP0INT(x, b) 

, b) - CHP0INT(x, a) 

b) - CHPOINTCx, b) - CHPOINTCx, 


a) 


I try applying G P Speck’s shading utility from 1998 with DERIVE 6.1: 


# 1 : 

# 2 : 


#3: 


# 4 : 


ISHADING1.MTH - first preload SHAD_UT . MTH| 


Concentric curve shading of a region bounded by parametric curve 
made up of portions of y = 1 + x A 2 and y = - x A 2 + 2x + 5 
XllCt) := i. ■ C H L E SS C t , -1) 

X12Ct) := t-CHCLOSEDCt, -1, 2) 

X13Ct) := C-t + 4)-CH0PENLFCt, 2, 5) 

X14Ct) := t-CHGRTCt, 5) 

YllCt) := i-CHLESSCt, -1) 

2 

Y12Ct) := Cl + t ) ■CHCLOSEDCt, -1, 2) 

2 

Y13Ct) := C- C-t + 4) + 2- C-t + 4) + 5) -CHOPENLFCt, 2, 5) 


YL4Ct) := i-CHGRTCt, 5) 



D-N-L#31 


G P Speck: Concentric Curve Shading 


P 5 


#5: 


Xl(t) := Xll(t) + X12(t) + X13(t) + X14(t) 
Yl(t) := Yll(t) + Y12(t) + Y13(t) + Y14(t) 

2 2 


#6: Ls(t) :: 

: 6-t - 

#7: APPEND 

" Ul ' 

s 

.. vl . 


2*x + 5, p:=0.5, q:=3] 

, APPEND( [ [Xl(S(t)) , YlCS(t))]], VECT0R([a* (Xl(S(t)) - p) + 


p, a.(Yl(S(t)) - q) + q], a, 0.1, 0.9, 0.1)) 


As you can see the plot is not really satisfying because of 
many gaps - which I cannot explain.. 

So I take the “alternative” IF-method - which caused prob- 
lems in various earlier versions. 

Unfortunately I cannot present an explanation for details like 
S(t) = 6t -1 , (-t + 4) replaces x, ... 



+ 2.x + 5 ] 


#1: Lul 1 + x , vl 

X2(t) := 

If -1 < t < 2 
#2: t 

If 2 < t < 5 
-t + 4 


Y2(t) := 

If -1 < t < 2 
#3: 1 + t A 2 

If 2 < t < 5 

- (-t + 4) A 2 + 2-C-t + 4) + 5 

#4: [S(t) := 6 ■ t - 1, u := X2(S(t)), v := Y2(S(t)), p := 0.5, q := S] 

ul 


#5 : APPEND 


vl 


, APPEND([[u, v] ] , VECTOR ([a- (u - p) + p, a-(v - q) + q], a, 0.1, 0.9, 0.1)) 


Original plot from 1998 (below) and 
DERIVE 6.10 plot (right). 0 < t < 1 for all plots 
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Comments from 1998: 

I tried to reproduce Mr Speck's graphs and found that the various DfW versions show different behav- 
iours: 

In version 4.0 the "Char"-functions don't work (Sorry, the highlighted expression cannot be plotted), 
and the "If" shows only the bottom half of the shading. 

Version 4.05 reacts on "Char" with "No finite real values!", "If" works. 

Version 4.10 similar to 4.05. Very strange behaviour with "Char" - after a few attempts mixes up the 
algebra window with the plot window, hangs up, 

I had no problems with DOS versions 3.10 to 4.x. Josef 

Al Rich's (SWHH) answer: Apparently the problem with the "Char"-functions is due to the huge ex- 
pression that is generated which in turn causes a memory overflow. We will investigate 

Since the IF( ) statement version of the shading program is faster and shorter, why publish the "Char" 
function version? 

Good question. Because it is an interesting example of "programming" with DERIVE. 


Concentric shading of the unit square given parametrically (IF-method): 

[GOO := F(t - 0.25), FOO := 0.5 + 0. 5 .C0S(4.?r.t)] U 


rX2(t) := 

If 0 < t < 0.25 
0 

If 0.25 < t < 0.5 

FOO 

If 0.5 < t < 0.75 
1 

If 0.75 < t < 1 

GOO 

[SCO := t, p := 0.5, q := 0.5] 

[u := X2CSO0), v := Y2CSO0)] 


Y2O0 := 

If 0 < t < 0.25 

GOO 

If 0.25 < t < 0.5 
1 

If 0.5 < t < 0.75 

FOO 

If 0.75 < t < 1 
0 


O.i. 

0.6 

0 > 

o.; 


0.2 0.4 0.6 0.8 


APPENDC[[u, v]], VECTORC[a.Cu - p) + p, a-Cv - q) + q], a, 0.1, 0.9, 0.1)) 


1 


Concentric shading of region made up of portions of bounding (curves FI, F2 and F3): 



APPEND([ [u , v] ] , VECTORCfa-Cu - p) + p, a-(v - q) + q], a, 0.1, 0.9, 0.1)) 



Concentric shading of a region made up of portions of bounding (curves FI , F2, F3 and F4): 

12 1 2 2 

FlOO := -x - ■ x - 2, F2(x) := 6 - x, F3(x) := x , F4Cx) - x - 2-x 

3 3 


[Y21(t) := F100, Y2200 := F2(S - t) , Y23(t) := F3(S - t) , Y24Ct) := F4CS - t)] 


rX2Ct) := 

If -2 < t < 4 
t 

If 4 < t < 10 
S - t 


Y2CO := 

If -2 < t < 4 
Y21(t) 

If 4 < t < 6 
Y22(t) 

If 6 < t < 9 
Y2300 

If 9 < t < 10 
Y24(t) 


[SCO := 12- 1 - 2, p := 1, q := 0] 

[u := X2CSO0), v := Y2CSO0)] 

APPENDC[[u, v]], VECTORC[a-Cu - p) + p, a-(v - q) + q], a, 0.1, 0.9, 0.1)) 

DERIVE 6.10 (below) vs DfW4 (right): 



5 



p 


8 


G P Speck: Concentric Curve Shading 


D-N-L#31 


Shading of one particular loop of [3 sin(t)+cos(0.5t), cos(5t)+sin(t)] with -n < t < +n. 
We fill the loop given by 0.14 < t < 0.84 and 2.07 < t < 2.72. 



Unfortunately I cannot reproduce G P Speck’s shaded figure because - shame on me - I 
cannot define the right IF-constructs. You can see the figure from 1998. 1 hope that George 
will explain the procedure more detailed and then we can try again. It might be a challenge 
for you, too: 



Memories of Gettyburg 


Carl Leinbach produced a CD-ROM containing the lectures and workshops of the 3rd International 
DERIVE and 77-92 Conference. You can purchase this inexpensive CD-ROM at Jerry Glynn's com- 
pany Math Ware. 


MathWare, P.O.Box 3025 Urbana, IL 61803, email: info@mathware.com 
Three keynotes can be downloaded from Carl's home page! 

As time went by this home page is not valid in 2013. 

One keynote is among the files collected in mth31.zip. 

You can download the Proceedings from: 

http://rfdz.ph-noe.ac.at/acdca/konferenzen/gettvsburg-1998.html 


andTI-92 

CONFE REN C E 
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Cyclone 98 


Now you are able to investigate graphs of surfaces which have simple equations which would 
be normally hard (or impossible) to graph like abs(x * y * z) = 1 or sin(2x) + sin(2y) = 
sin(2z). You can graph 3D inequalities and you also can generate lists of surfaces. The objects 
can be animated. Cyclone 98 runs on Win95 or Windows NT. 


I googled for Cyclone 98 - and what did I find? There were al lot of advertisements for two 
products: 


A paintball equipment - 


and a model of a famous race car. 




But there is no soft ware!! Go for http://www.dpgraph.com/ 


Philip Yorke asked me to print the following paragraph: 

Chartwell-Yorke have added full details of a wide range of English language books using DERIVE to 
their website. You will find full descriptions, cover photos and contents lists at 
http://ChartwellYorke.com/Booksl.html. There are a great many new books which you may not be fa- 
miliar with on the site! They would like to add your independent comments about the books to your 
site, so you can help each other to choose which suits your needs. If you are family with any of the 
books covered, please email your comments on their usefulness to Philip.Yorke@ChartwellYorke.com 
and he will add your feedback to the site. 
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From ' 

"nested IF expressions" to a "Macro for Solving Equations" 

From: 

HeinzRainer.Geyer@t-online.de 

31-01-1998 

Subject: 

nested IF expressions 



Hi all, what is wrong with the following definition: 

TRANS F (equation , operator, value ) : =IF (operator = M +" , equation+value, 

IF (operator = equation-value, "?") ) 

Test TRANSF ( 1 Ox + 6=0 , " - " , 6 ) returns the unsimplified expression, 

with TRANSF ( 1 Ox + 6=0 6 ) you get what one is expecting: 1 Ox = 6 

I want to define a function to automate equivalence transformations of equations, of course with more 
than 

I use DfW 4.06 German and DfW 4.0 English. 

Thanks for hints. Helpful would also be an other form than the nested if structure. 

Rainer Geyer 


From: math@rhombus.be 31-01-1998 

Dear Heinz, 

I can't solve your nested IF's problem, but: 

please see the attached LINEQ.MTH file which automates and illustrates the consecutive equivalence 
transformations on linear equations. 

Best regards, Jan Vermeylen 

LIN_EG.MTH - Jan Vermeylen, version 4: 

#1: EQUIVO(eq) := EXPAND(eq) 

EQUIVl(eq) := EXPAND(eq) - (TERMS(LHS(EXPAND(eq)))) - (TERMS(RHS(EXPAND(eq)))) 

# 2 : 2 1 

EXPAND(eq) 

EQUIV2(eq) := 

#3: (FACT0RS(LHS(EXPAND(eq)))) 

2 


#4: 


#5: 


example 


x x - 2 

2 ■ (x - 1) + = 4 -x - 9 + 

3 3 


r x x — 2 

#6: EQUIV0 2-(x - 1) + = 4-x - 9 + 

L 3 3 j 


# 7 : EQUIV1 


7-x 


13 -x 29 


2 = 


7-x 
V 3 
23 > 


13 «x 


2 = 


2-x = 


#S: EQUIV2 


23 


2-x = 



29 ^ 


3 } 
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DERIVE 6.1 0's answers: 


#9 : EQUIVO 

13- x 

#10: x = 


x x - 2 ^ 

2-(x - 1) + = 4 -x - 9 + 

3 3 j 


23 


7 7 

Expand does much more than only expanding both sides!! 

#11 : EQUIV1 


' 13- x 23 ^ 


f 

v — 

28 ^ 

1 

■*vl 

1 

■*vl 


jS — 

6 j 


From: Spwelke@aol.com 01-02-1998 

Hallo Mr. Geyer, hallo all DERIVE users, 

Your problem is that DERIVE does not evaluate string expressions like "Hallo" = "hallo" as false, but 
DERIVE evaluates "Hallo" = "Hallo" as true: 

#1: IF ( "Hallo" = "Hallo", 1, 2) = 1 

#2: IF ("Hallo" = "hallo", 1, 2) = IF("Hallo" = "hallo", 1, 2) 

doesn't evaluate as 2 as one would expect. A simple way around your problem is to assign numerical 
values to the operations in mind, for example: 


DERIVE 6 works as expected: #2: if ("Hallo" = "hallo", l, 2) =2 


But we meet again the same problem as 
before, DERIVE does too much in this 
case. 

DERIVE immediately solves the equa- 
tion for the variable and does not con- 
sider the given task. 


transf(eq, op, val) 
If op = M + M 
eq + val 
If op = M - M 
eq - val 
eq 


transf (4 -x + 3=8, 



transf (4 -x +3=8, +, 300) 





There is a much simpler way to deal with #6: 
your particular problem. If I understood your 
intentions well, the problem is to handle #7: 
equations as we all have learned at school. 

# 8 : 

Please see an example: 

#9: 

We can see that arithmetic operations and 
functions apply to both sides of an equation 
separately. #1 “ 

# 11 : 


2 — x 

2 + 3 = 9 

2 — x 2 — x 

(2 + 8 = 9) - 8 = (2 =6) 


2 - x 

LOG (2 = 6) = ((2 - x) ■ LN(2) = LN(6)) 


(2 - x) . LN(2) = LN(6) 

f LN(6) ‘ 

25 n . H - 

LN(2) 

LN(6) 

f 

- 2 = -: 
s 

■c-i) = 

/ 

^ LN(2) . 

LN(6) h 

i. - -i 

a LN(2) j 

( LN(6) 

n . H - 25 

LN(2) j 

f LN(6) ■ 

TH - 25 

^ LN(2) 

^ LN(2) . 
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From: GNSStepan@aol.com 02-02-1998 

Hallo Mr Welke, hallo all DERIVE users, 

your simple way handling equations was very interesting to me. At school it is important that simplifi- 
cation occurs step by step. Executing your example with DfW 4.05 equation #10 simplifies to 


LN (3) 


1 instead of 


LN (6) 


LN ( 2 ) LN ( 2 ) 

What can I do to avoid this additional simplification step? Thanks for a tip. Georg Stepan 


(A short additional comment: You might give the problem a new quality: ask your students for 
DERIVE's hidden calculation, Josef) 

From: Spwelke@aol.com 02-02-1998 

Dear Mr. Geyer, dear Mr. Stepan, dear DERIVE users, 

(1) At first I respond to Mr. Geyer "s request concerning a CASE-Function: 

A sequence of nested IF-expressions calls for a CASE-function, which allows the user to select 
between more than two alternatives. A simple way to achieve this functionality is the following 
CASE-function: 


# 1 : 


#2 


CASE(a_, v_) := (SELECT (a_ = k_ , k_, v_)) 

1 


1,2 


test 


3 4 

4 5 

5 t 
7 z 

#3: CASE (5 , test) = t 

#4: [plus := 1, minus := 2, times := 3, dividecLby := 4] 




r 

" 1 

eq + val 





2 

eq - val 


#5: 

TRANSECeq, op, val) := CASE 

op, 

3 

eq ■ val 





4 

eq 

val 

j 


DfW 4: 


#7: TRANSE 


#S: TRANSE 


2-x 
v 3 
{ 2-x 
k 3 

f 

2-x = 


- 7 


11 


plus, 7 


2-x 


79 


79 

, times, 3 

11 j 


v 3 11; 

237 L 

2-x = 

V 11 J 


237 


11 


divided_by, 2 


> 


" 237 ^ 

2 

- 

'A - 

j 


L 22 J 


#9: TRANSE 
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but DERIVE 6 solves - unasked - the equation: 


# 10 : 


TRANSE 


2-x 


3 


- 7 = 


2 


11 


> 


r 237 ^ 

7 

- 

x = 

j 


L 22 J 


Another use of the CASE-function is to assign the value true or false to the first argument of 
the CASE-function. Then the first element of each line v_ must be a logical expression, which 
can be evaluated to true or false. Please see the next example: 




f 

< 

1 

too small 


#11: 

FOO := CASE 

true , 

> 

3 

too big 





_ X_ > 1 A 

x_ < 3 

ok 

j 


#12: [F(-3), F(e), F(2-rr)] = [too small, ok, too big] 

Works also in the recent versions, Josef 


(2) Now I respond to Mr. Stepan's question 

There is no way to manipulate the evaluation of expressions in DERIVE. In MATHEMATICA , for 
example, there is an extensive set of commands to control, prevent and to alter the way 
MATHEMATICA evaluates expressions. I cannot see any simple trick to prevent the evaluation of 
LOG(6)/LOG(2) as LOG(3)/LOG(2) - 1. Sorry! Maybe someone else in the DERIVE community 
has a trick? 

So long for this time, Stefan Welke 

The next part of this long way to the final Macro is of historic interest. I reprint the evolution 
process of this equation solving macro designed and improved by Johann Wiesenbauer. 

As the DERIVE versions changed his MACRO has to be changed, too. Commands which 
worked in earlier versions did not work any longer... You must consider that in 1998 we could 
not write complex programs as we can do now. There were “only” mor or less complicated 
functions which called other functions which called ... 

The final version which works even with DERIVE 6.10 was published in DNL#53 (March 
2004). I don’t intend to let you waiting until DNL#53 will be revised ©. I include the latest ver- 
sion in this DNL#31 rev. But you may try providing a similar MACRO for Tl-Nspire, Josef 

From: J.Wiesenbauer@tuwien.ac.at 03-02-1 998 

Subject: Once again: Nested If expressions 

Hi all, 

Yes it is true, the idea of solving an equation step by step is great, but I have not yet seen a really satis- 
fying solution. What I think is the best approach makes use of the following tiny auxiliary function: 

#1: DO (u) := (eq := u) 

In the following examples you can see it at work: 
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Example 1 (Jan Vermeylen's equation) 


#2: 

eq 

:= 2- 

• (x - 

- i) + 

X 

— ZL • v 

X 

o + 

- 

2 

3 

4 A 


3 


#3: 

DO (EXPAND (eq) ) = 

r 7< 

• X 

D 

13 • x 


29 

L 3 

3 


3 

#4: 

DO 

[eq - 

13- 

>x 1 

_L D - 

■ [- - 

23 

1 


3 -1 

3 

J 


#5: 

DO 

r eq 1 

- 1 

r„ _ 23 1 





L -2 J 

1 


6 J 





#6: 

"or : " 








#7: 

eq 

:= 2- 

• (x - 

- 1) + 

X 

— ZL • v 

X 

Q _i_ 

- 

2 

3 


Z> \ 

3 


#8: 

DO (EXPAND (eq) ) = 

r 7 ' 

• X 

D 

13 * x 


29 

L 3 

3 


3 

#9: 

DO (eq- 3) 


(7 • x - 

6 = 

= 13 • x - 

29) 



#10: 

DO(eq - 

13 • x + 6 ) 


( - 6 • x = 

-23) 



#11: 

DO 

r eq 1 

- 1 

r 23 i 





L -6 J 

■ 1 

L 

6 J 






Assets: (1) My method is not restricted to linear equations (cf. Ex. 2) 

(2) There is no need to input equarions (except for #2)! 

Remark: In contrast to Jan's method the pupil has to find the appropriate transformations by himself. 

In my eyes this is actually an advantage (after all, we know that DERIVE can do it!), but 
you may view this as a draw-back 

Example 2 (Stefan Welke's equation) 


#12: 

eq : = 2 

+ 3 

#13: 

DO (eq - 3) 

2 

= (2 

#14: 

DO (LOG (eq. 

2) ) = 

#15: 

DO (eq - 2) 

- b 

#16: 

r e T 1 

D °[ ] = 

[*- 


= 6 ) 


LN (3) 


+ 1 


LN (2) 

-4 


LN (3) 

LN (2) 

LN(3) -I 
LN ( 2 ) -I 

Remark: At first glance this is very close to Stefan's solution, but unlike him I don't have to care 
about equation numbers and, above all, what looks nice in the input line at first turns immediately into 
a very ugly term after pressing ENTER! What you see (in his e-mail) is not what you get (on the 
screen)! 

(In his e-mail Stefan used the expression numbers as an abbreviation. Johann closed his message 
with another function to give a concluding report containing all steps performed. I'll skip that part be- 
cause he overtrumped himself in the message which arrived 10 days later. Josef) 

Don’t try now, because the “tiny auxiliary” function DO(u) does not work!! 
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From: J.Wiesenbauer@tuwien.ac.at 13-02-1998 

Subject: A macro for solving equations 

In the attached file you can see for the first time what one might a call a DERIVE-macro\ I have used 
it to record the step-by-step solution of an equation, but advanced programmers should have no diffi- 
culty in exploiting the underlying ideas for other purposes! (Good luck!) 

Cheers, Johann 

#1: "A macro for solving equations ((c) Johann Wiesenbauer, Vienna) " 

#2: "The following six functions must be loaded:" 

#3: NEWO (u) := (o:=u) 

#4: NEWEQ(u) :=(eq:=u) 

#5: DO (u) :=NEWEQ (LIM (u, eq) ) +0*SUM (NEWO (APPEND (o, [u] ) ) ) 

#6 : UNDO (u) :=NEWEQ (LIM (ITERATE (u,@,@,-l) ,0,eq) ) +0* SUM (NEWO 

( DELETE_ELEMENT (o, DIMENSION ( o ) ) ) ) 

#7: RECORD (u) : = (eq: = (e : =u) )+SUM(o:=[] ) 

#8: play:=DELETE_ELEMENT (ITERATES ( [LIM (a_ SUB l,0,e_),IF 

(DIMENSION (o_) =0, o_, [o_ SUB 1 ] ) , DELETE_ELEMENT (o_) ] , [e_, a_, o_] , 
[e, [o SUB 1] , DELETE_ELEMENT (o) ], DIMENSION (o) )', 3) ' 

#9: "Study the following examples carefully to see how it works!" 

2 • x + 1 

#10: RECORD (3 - 7 = 20) User 

2 • x + 1 

#11: 3 - 7 = 20 S imp (#10) 

#12: DO ( @ + 7) User 

2 • x + 1 

#13: 3 =27 S imp (#12) 

#14: DO (LOG ( @ , 3)) 
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I skip the 2 nd example, because this version does not work in recent versions. Josef 

From: J.Wiesenbauer@tuwien.ac.at 14-02-1998 

Subject: Troubleshooting as regards my macro 

Since I got e-mails of users who have troubles with my macro two comments might come in handy: 


• All evocations of my functions must be simplified subsequently, in particular RECORD ( . . .). 
Sorry that I was mistaken in taking that for granted. (After all, you can see those simplifications in 
the annotations of my example! ! !) 

• If you still have troubles, this might be due to an older version of DfW. I can only guarantee that it 
works in DfW 4.08!) Just in case you have an older English version DfW 4.0x, there is an easy 
remedy: You should use a free (!) update on the DERIVE- homepage! ! ! 

• Cheers, Johann 

(I tried to use Johann's macro with DOS versions 3.x. It works if you replace all the @s in 
do, undo and play by any unused variable, e.g. take g_ . I wish much luck, much fun and 
much successworking with this macro. Josef) 


This is Johann’s the final MACRO - version (together with small additions. Johann wrote: 

“On the other hand, I have added a new command sqr(@) to compute square roots. In a 
similar way, other useful commands (e.g. for expanding and factoring polynomials) could be 
added, if there is a need for it." 

Examples: do(@,”e”), do(@,”f ) and do(@,”r”) for expanding, factoring and drawing the root. 


RECORDfiO := 

Prog 

0 := □ 

@1 := LHSfiO 
#1: @2 := RHSfiO 

1 := [@1] 
r := [@2] 
u 

sqrfu, O := 

Prog 

If NUMBER?(iO 
#2: RETURN ± ^u 

u_ := P0LY_GCD(u, 3(u, x)) 
If NUMBER? Cu/u_ A 2) 
•/Cu/u_ a 2) ■ u_ 

UND0(l0 := 

Prog 

0 := REST Co ) 

1 := REST OJ 

#4: r := REST CD 

@1 := FIRST CD 
@2 := FIRST CD 


doCu, t := 0 , t_) := 

Prog 

t_ := CNAME_TO_CODESCtDll 
If t_ = 101 
Prog 

@1 := EXPAN DC@1 , '’Radical") 

@2 := EXPANDC@2, "Radical") 
o := ADJOINC["e>!pandC@)" ] , o) 
If t_ = 102 
Prog 

@1 := FACTOR C@l, "Radical") 

@2 := FACTOR C@2, "Radical") 
o := ADJOINC["factorC@)" ] , o) 
#3: If t_ = 114 

Prog 

@1 := sqrC@l) 

@2 := sqrC@2) 

o := ADJ0INC['W@"], o) 

If t_ = 48 
Prog 

@1 := SUBSTCu, @1) 

@2 := SUBSTCu, @2) 

o := ADJOINCtu], o) 

1 := ADJOIN C@1 , D 
r := ADJOIN C@2 , r) 


# 5 : 


play := [REVERSE(l) , VECT0R(=, k_, DIM(l)), REVERSE(r), APPEND(REVERSE(o) , [])]’ 
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A few examples: 


#6: RECORD^ (3 


3 - / 


#8: do(@ + 


#9: do (LOG (@ , 3)) = 


r @ 

#10: do 


#25: REC0RD( (x - 4-x + 7 = 10)) 


#26: x - 4-x + 7 = 10 


#11: Oops, sorry! 


#12: UNDOO 


#13: do(@ - 1) = 

#14: dof— ] = [ 


3 - / 


[@ + 


#15: play = 2-x + 1 = 3 [0 - 1] 


#27: do(@ 


#28: do( JiS) = 


- 4-x + 4 = 7 


This is the standard / let's try the parameter M r M 


#29: do(@, r) = 


#30: do(@ + 2) = 


x - 2 = ± yV 


A more complex equation: 


= + J7 + 


f'f b + x 4 ■ b ■ (3 ■ b - x) b - x Y| 

RECORD = 

#16: b - x 2 2 b+x 

^ L b - 2-b-x + x )j 

b+x 4-b- (3-b - x) b - x 


#17: b — x 


b - 2-b-x + x 


#18: do(@-(b - x) ■ (b + x)) = 


(x + b)-(x - 4-b-x + 11 -b ) = (b - x) 


#19: do (@ , e) 


? ? ^ 


-b-x - 7 ■ b ■ >: - 11- b = - x + 3-b-x - 3-b ■ >: + b 


#20: do(@ + x — 3-b-x + 7-b -x) = 3 23 

_ - 11 -b = 4-b-x+b 
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This is the last lecture of the 1997 International DERIVE and TI-92 symposium held in Saro, Sweden, 
to be published in the DNL. 


An Algorithm for Extracting Logic Propositions from Numerical Data 
and its Implementation with DERIVE™. 


Terence Etchells 

School of Computing and Mathematical Sciences 
Liverpool John Moores University, UK. 

Abstract 

This article is based on the work of Hiroshi Tsukimoto of the Systems and Software Engineering 
Laboratory, Toshiba Corporation [1]. The aim of this work is to review the work of Tsukimoto 
and implement his algorithms in Derive for Windows. 

If a regression analysis has been performed on some data set, then the resulting expression is an 
algebraic one involving variables x, y, z, ... etc. depending on the number of variables. For exam- 
ple, imagine that we are given a set of data of three varaibles x, y and z- A linear regression plane 
of the form z = ax + by + c can be readily obtained. What does this regression line tell us about 
the underlying structure of the data? What we require is a Logic proposition that describes the 
hidden structure that generates the data. For example such a proposition could be z = (x a — i y ) , 

which would be interpreted as the property z is true if the property x is true and the property y is 
not true (i.e. false). The essence of this paper is to describe and implement an algorithm for ex- 
tracting such propositions from given regression hyper planes or data. 

Backgound 
Logic Vectors 

Propositions in Boolean logic can be represented as vectors, theses vectors are called Logic Vectors. 
We will spend some time looking at the construction of Logic vectors and their building blocks at- 
oms. As it was thought for some time that atoms were the building blocks of matter, atoms in the 
Boolean sense are the building blocks of Boolean functions. 

Let us take the most simple Boolean function x. The function x can take only two values 0 or 
1 (i.e. elements from the set {0,1}) as it is Boolean, and also consider its negation — iX which can also 
only take values {0,1}. The truth table of x and — iX is 


X 

— iX 

1 

0 

0 

l 


Figure 1 

Truth table for -ix 


which tells us some thing very important, that all Boolean functions of a single variable will either be 
a contradiction (e.g. x a — iX ), a tautology ( e.g. x v — ix ) , x or —oc . This is because a truth table of 
a Boolean function of one variable can only contain the values [1,1], [1,0], [0,1] or [0,0]. Here we have 
written the vertical column of the truth table as a horizontal vector, the idea of the truth values of a 
function being represented as a vector is paramount to the ensuing work. To illustrate this further we 
will study the truth tables of some Boolean functions of a single variable, namely 
— i((x a — ix) v — i(x ax)); — iX v — i(x ax); x a x v x . 

As DERIVE is used extensively in this work we will use it from the outset to produce truth tables. 
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TRUTH_TABLE(x , ((x a x) v (x a x)), x v (x a x), (x a x) v x) 

X C(x A x) V (x A x)) -iXV-i(xAx) (xAx)vx' 

true true false true 

_ false false true false 

As it is desirable to have values of 0 or 1 instead of the variables true and false, the functions 
TEST(arg) and TRUTH (tt) are introduced. 

TEST (arg) := 

If arg = true 
1 
0 

arg 

TRUTH(tt) := APPENDCT tt 1, VECTOR (VECTOR (TEST (tt ), m, 1, DIMENSI0N(tt 1 )) , r, 2, DIMENSION(tt))) 

L lJ r.m 

X -.((XA--X)V-.(XAX)) -.XV-i(XAX) (XAX)vx' 

TRUTH true true false true 

_ false false true false 

x x -.x x 
110 1 
_ 0 0 1 0 _ 

As we can see the application of the TRUTH () function simplifies the Boolean expressions on the top 
row of the truth table to their minimised (reduced) form. DERIVE ’s ability to reduce Boolean expres- 
sions to a minimum form automatically is utilised in later work. So any Boolean function of a single 
variable can be written in terms of the ‘atoms’ x and —oc . Indeed, they may reduce to a tautology 
[1,1] or a contradiction [0,0], but these are not atoms because an atom must be orthogonal to all other 
atoms and be of unit size. We should think of these ‘atoms’ as orthogonal unit vectors in a 
2 dimensional Euclidean space, and as we do not think of matter comprising of nothing, or indeed 
everything, the Logic vectors [0,0] and [1,1] are not atoms. 

We now extend these ideas to Boolean functions of more than one variable. The atoms of single vari- 
able functions arose very naturally from the truth tables, so the truth tables of Boolean functions of x 
andy should help the development of the 2 variable atoms. 

" x 1 1 0 0 " 

TRUTH (TRUTH_TABLE(x, y)) 1 = 

_ y 1 0 1 0 _ 

As it can be seen from the truth table above, for Boolean functions of two variables, the vector that 

represent their truth values have dimension 2 2 = 4 . We can see that the truth vector of x = [1,1, 0,0] 
and that ofy = [1,0, 1,0], so neither x nor y are atoms in a 2 variable Boolean expression as they are not 
of unit length. Hidden in the truth table above is the structure of the atoms! Reading a 1 as a proposi- 
tion and a 0 as its negation the truth table above can be read as; 

[x a y , x a —i y, — iX a y , — oc a — y] 

We will now verify that x a y, x a — y, —oc a y and —oc a — y are indeed the atoms of Boolean func- 
tions of two variables x andy. 
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The next few lines of DERIVE commands show the linearisation 

( 3 3 2 2 \ 

x y +x y +xy + x + y + \) = 3 xy + x + y + 1 in action. 


3 3 2 2 

REMAINDER (x -y + x -y + x -y + x + y + 1 , x - (1 - x) , x) 

3 2 

x ■ (y +y +y+l)+y+l 
3 2 

REMAINDER (x ■ (y +y + y + 1) + y + 1, y.(l - y), y) 
x-(3-y + 1) + y + 1 

The same expression in a slightly different form. We will now automate this process with user defined 
DERIVE functions. An extremely useful DERIVE function is VARIABLES (u), which returns a vector 
of all the free variables in the expressions u, e.g. VARIABLES(3xy+2x+3z) returns [x,y,z] when 
evaluated. The next two functions automate the linearisation procedure. 

t_AUX(wlw, vlv) := (ITERATE(T rlr + 1, REMAINDER(ulu , vlv • (1 - vlv ) , vlv )], 

L rlr rlr rlr J 


[rlr, ulu], [1, wlw] , DIMENSION(vlv))) 

2 


t (wlw) := t_AUX(wlw, VARIABLES(wlw)) 

Programming notes: 

rlr, wlw, vlv are obscure variables to avoid clashes with functions that may be defined later and call 
this particular function. It is not wise practice to use the same variables in a functions that call one an- 
other. The V SUB r function extracts the r th element of the vector V. The use of auxiliary functions 
avoids the repetition of calculations, by passing evaluated results into other functions. 

The x function in action: 

3 3 2 2 

t(x -y + x -y + x-y + x + y + 1) = x-(3-y + 1) + y + 1 

We now have a function that will automatically linearise a polynomial of any order and any number of 
variables. Functions of this type, with two or more variables are called multi linear functions. 

The x function can be used to simplify a sequence of a operations on many variables, for example the 
Boolean expression ((x axaj/)a(xaj/a z)) a(xaz) converts to the elementary algebraic ex- 

4 2 2 4 2 2 

pression (( xxy )(xyz))( xz ) = x y z and t(x y z ) = x y z , which converts back to x a y az . 
More formally we have 


x a y <=> z(xy) . 










p 24 T.Etchells: Extracting Logic Propositions from Numerical Data D-N-L#31 


n 

<f>t = Y\ e { x j’ BD( j~ l ’ n + W)) (^ = I ••• , 2 ", 7 = 1 , ••• , «) 

j = i 

we will elaborate on how this definition works. 

In the example above, we used n = 3 variables and this generates the orthonormal expressions 


orthonormal 


expressions 

p values 

x l x 2 x 3 

000 

* 1 * 2 ( 1 - * 3 ) 

001 

*1 (1 *2 )*3 

010 

*1 (! *2 )(l *3 ) 

Oil 


we can see that the switch to decide whether p = 0 or 1 forms a binary sequence starting at 0 and finish- 
ing at 2 n - 1 . 

So that this definition may be coded we need a function that will convert a decimal integer in to its bi- 
nary equivalent. To make use of the binary digits we will need to produce a vector of the binary dig- 
its, e.g. convert 5 to the vector [1,0,1]. The function BINARY_AUX(m2m, n 2 reproduces such a vector 
of the first m2m binary digits of the decimal integer n2n. This function relies on the division by 2 algo- 
rithm that is taught (or used to be) in most high school mathematics curricula. 


BINARY_AUX(m2m , n2n) ;= R EVE RSE_VECTOR ITERATES 


p2p - M0D(p2p, 2) 
2 


T n2n - M0D(n2n, 2) “I h 

[p2p, r], , M0D(n2n, 2) , m2m - 1 ' 

2 J J 

We now define the functions 

EEl(xlx, pip) := 

If pip = 0 
xlx 

1 - xlx 
1 - xlx 


, M0D(p2p , 2) 


<J)(vars, ill, nln) := FI EEl(vars , (BINARY_AUX(nln , ill)) ) 

j=l j j 

to automate the above definitions. 

To see these functions in action, we will find the orthonormal expressions for the variables [x,jy,z] . 
VECTOR (cp([x, y, z], i, 3), i, 0, 7) 


[x-y-z, x-y- (1 - z), x-z-(l - y), x-(y - l)-(z - 1), y-z-(l - x), 
y.(x - l).(z - 1), z-(x - l).(y - 1), (1 - z).(x - l).(y - 1)] 
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To automate the process of producing a vector of the orthonormal expressions we define the functions: 

num 

ORTHONORM_AUX(exps , vars, num) := VECTOR(c|)(vars , i, num), i, 0, 2 - 1) 

ORTHONORM(exps) :z 0RTH0N0RM_AUX(exps , VARIABLES(exps) , 

DIMENSION (VARIABLES(exps))) 

Our task is to convert an expression, say xy — x + l, into its orthonormal form. We can do this by 
writing, for example, 

xy - x + 1 = axy + bx( 1 -y) + c( 1 - x)y + d( 1 - x)(l - y) 

Clearly, we can find the coefficient a by substituting x = 1 and y = 1, (1,1) into the above expression. 
We get 

1 • 1 — 1 + 1 — fl*l*l + Z?*l*0 + c*0*l + t/*0*0 ci — \ 

repeating with (1,0) ,(0,1) and (0,0) will give b , c and d respectively. Notice that we are substituting 
the coordinates of a square. In general, for functions of n variables we substitute the coordinates of an 
n dimensional hyper-cube into the function to find the elements of the logic vector. 

The function HYPER_CUBE(n3n) produces a matrix of coordinates of the n3n unit hyper-cube. 

n3n 

HYPER_CUBE(n3n) :z VECT0R(BINARY_AUX(n3n , w3w), w3w, 0, 2 - 1) 

The following functions calcuates the Logic vector of a Boolean function in T logic form. 

L0GIC_VECT0R_AUX(u2u , vars, vect) :z VECT0R( lim u2u, r2r, 

vars-^vect |r2r 

DIMENSION(vect) , 1, -1) 

L0GIC_VECT0R_AUX2(u2u , vars) :z L0GIC_VECT0R_AUX(u2u , vars, 
HYPER_CUBE(DIMENSION(vars))) 

L0GIC_VECT0R(u2u) := LOGIC_VECTOR_AUX2(u2u , VARIABLES(u2u)) 

Using this function on xy - x + 1 we have; 

L0GIC_VECT0R(x.y - x + 1) = [1, 0, 1, 1] 

Logic Propositions from Regression Functions 

Until the work of Tsukimoto and Morita there have not been any methods for extracting logical pro- 
positions from regression planes. For example, the function 

z = 0.71x-0.99y + 0.25 

is obtained from a least squares fit of a set of data points. What is the underlying Boolean structure of 
the data. Looking at the equation, if we consider in the region [0,1] near 0 is not true and near 1 is 
true. If we apply our method for finding Logic vectors to this regression line we will get a vector 
whose elements are in the interval [0,1] instead of the set {0,1}. 
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NEGATE_MAT_AUX_l(www) := NEGATE_MAT_AUX(www, DIMENSION(www)) 

NEGATE_MAT(ff) NEGATE_MAT_AUX_l(VARIABLES(ff)) 

Produces an n x 2 matrix with the n variables of f f in the first column and their Boolean negation in 
the second column. 


VECT_AND_AUX(vars , nnn) 

:= (rTERATES(Trrr + 1, vvv a vars "1, 


L rrr + lj 

[rrr, vvv], |“ 1 , vars 

, nnn - 1)) 

nnn, 2 

VECT_AND(vect) := VECT_AND_AUX(vect , DIMENSION(vect)) 

Converts a vector [x, y, z, .. 

. . ] into the logical expression x a y a z a ••• 


VECT_OR_AUX(vars, nnn) :z (rTERATES(T rrr + 1, vvv v vars 

rrr + 1 


[rrr, vvv], |~1, vars j, nnn - 1)) 

nnn, 2 

VECT_OR(vect) :z VECT_OR_AUX(vect , DIMENSION(vect)) 

Converts a vector [x, y, z , . . . . ] into the logical expression x v y v z v ••• 

ZERO_FALSE(vect) := VECTOR (IF (vect = 0, false, vect , vect ), r, 1, 

r r r 

DIMENSION (vect)) 

Converts any 0’ in a vector to the statement false. 

r 

PR0P0SITI0N_AND_AUX(ff , nmat) := VECTOR VECT_AND(nmat ), r_, 1, 

< r_ 

DIMENSION(VARIABLES(ff ))) 

2 

y 

PR0P0SITI0N_AND(ff) := PROPOSITION_AND_AUX(ff , NEGATE_MAT ( f f ) ) 

Creates a vector of all atoms for ff, e.g if ff=xy+x then this function produces the vector 
[x a yjc a— iy ,— cca y,— a a— i_y] 

ELEMENT_PROD(vectl, vect2) := VECTOR(vectl -vect2 , r_, 1, 

r r 

DIMENSION(vectl)) 
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Produces an element by element vector product (Logic And). 


EXTRACT_P ROP ( f f ) := VECT_OR(ZERO_FALSE(ELEMENT_PROD(MAX_LIKELY_HOOD(ff ) , 
PROPOSITION_AND(ff)))) 


This final function takes any multilinear function and produces a pseudo-maximum likelyhood Boo- 
lean function that best describes the underlying Boolean structure. 


Experiment 

These functions will generate ‘noisy’ data that has an underlying Boolean structure of x a —>y a z. 

F3(x) := , G3(y) := , H3(z) := 

If x < 0.5 If y < 0.6 If z < 0.4 
0 0 0 

1 1 1 

PP(tlt) := 

If tit = 0 
RANDOM(l)/2 
RANDOM(l)/2 + 1/2 
RANDOM(l)/2 + 1/2 

APPEND(APPEND(VECTOR(VECTOR(VECTOR([x , y, z, PP(F3(x) v (A G3(y) a 
H3(z)))], x, 0, 1, 0.25), y, 0, 1, 0.25), z, 0, 1, 0.25))) 


Using DERIVE’s FIT() function we can find a multilinear regression curve of the form 
w = axyz + bxy + cxz + dyz + ex + f y + g z + h . In this experiment it turns out to be 

0.932032-x-y-z - 0.0892345-x-y - 0.817753-x-z + 0.629476-x - 0.807328-y-z + 

+ 0. 103246 -y + 0. 724527 -z + 0.188913 

We now use our function to exctract the underlying Boolean proposition 
EXTRACT_PROP( 0.932032-x-y-z - 0.0892345-x-y - 0.817753-x-z + 0.629476-x - 
0.807328-y-z + 0.103246-y + 0.724527-z + 0.188913) 


x v y a z) 


Success!! 


References 

[1] The Discovery of Logical Propositions from Noisy Data, Hiroshi Tsukimoto, AAAI Work- 
shop on Knowledge Discovery in Databases, pp 205-216. 


Comment for DERIVE 6.10: The latest version of DERIVE confirms all intermediate results. 
An exception is the final outcome because of the leading zeros in the logic expression: 

0v0v0vxv(->yAz) 

x v (-> y a z) v 0 


Hi Terence, can you explain and adjust the respective function(s)? 
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Dear Editor, 

We are Ukrainian followers of Information Technologies (in particular DERIVE and CABRI, unfortu- 
nately not TI-92 due to financial circumstances) in Mathematics education for many years. We are 
now subscribers of your remarkable journal. In this matter we send you our article which maybe inter- 
esting for publishing in your journal. 

With best regards to you and your journal, 
sincerely yours, 

Dr. Rakov S.A. and Dr. Gorokh V.P. 

Information Technologies in Geometry 

(An example of the generalization of a well-known problem about squares) 

Rakov S.A . & Gorokh V.P. 

Kharkov State Pedagogical University named after G.S.Skovorada 
tel/ fax : (+380) (0572) 676-604 
e.mail: Rakov@pu.ac. Kharkov. UA 


Introduction 

Information technologies have changed all kinds of human activities. Mathematics is not the exception 
- it becomes technologically dependent. The most important changes take place in the process of doing 
mathematics - discovering new facts and their proof. Special mathematical packages offer the user a 
suitable environment for arranging computer experiments in the problem field with the aim of finding 
mathematical regularities on the first step of exploration and then support the process of proof with the 
powerful opportunities of computer algebra. No doubt that the future of mathematics is a symbiosis of 
man and computer. Using mathematical packages becomes the inalienable component of mathematical 
culture. 

Innovative trends in mathematical education are lying in the framework of constructive approach - in- 
volving students in the process of constructing their own mathematical system which consists of 
mathematical knowledge and mathematical beliefs. One of the most effective ways of the realization 
of a constructive approach is the method of learnig explorations when students explore open-ended 
problems on their own. Solving open-ended problems can be regarded as the model of the professional 
mathematical work. Therefore it is naturally to use information technologies in mathematical educa- 
tion just in the same manner: computer experiments as the source of powerful ideas, computer algebra 
as a tool of deductive method. Using information technologies in arranging learning explorations and 
carrying out the proofs cannot only do this work more effective but acquaint students with modem 
technologies of mathematics. 

The authors have had the pleasure of playing with Cabri-Geometre in one well-known problem [1] and 
as a result have conjectured its generalization which then they have proved by the help of DERIVE. 
They were highly impressed with these computer games and decided to share their experience with 
their students. Thereby the 2-week computer practice with undergraduate students at the Mathematical 
Teachers Department of Kharkov State Pedagogical University in 1997 was organized in the following 
form: 1st week - acquainting with packages Cabri-Geometre and DERIVE , 2nd week playing with the 
above mentioned problem with the aim of its generalization. The students have coped with the task 
posed successfully without any help. The current article outlines activities which were used in the 
process of generalization of this problem (it is rather common for authors and students). 
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Problem 

The problem area is: 



Two sqares A X QC X D X and A 2 B 2 C 2 Q with centers O x and 0 2 and common vertex Q with identical orien- 
tation are given. Let E and F the midpoints of the segments A X A 2 and CiC 2 . Then the quadrangle 
0\E0 2 F is a square as well. 

How can the problem mentioned above be generalized? 


Computer experiments in solving the problem 
Experiment 1 

At first we construct a computer model of a given problem in the Cabri-Geometre environment. For 
that purpose we first construct a macro-construction quad which draws a square using two given 
points at the endpoint of the diagonal. 

(The authors presented Cabri constructions. I reproduced the constructions with the Tl-92's Geometry 
Tool which is a Cabri implementation on the TI-92. You can find the .92a files on the Diskette of the 
Year. Josef) 

Produce a computer Model 1 of the problem: 

(1) Construct a macro quad which draws a square on two given points at the ends of the diagonal. 

(2) Define three points D x , Q and B 2 and draw the squares A X QC X D X and A 2 B 2 C 2 Q using pairs of 
points D\, Q and Q , B 2 as endpoints of the diagonals to apply the macro quad. 

(3) Draw the segments A X A 2 , QB 2 , CiC 2 and D X Q. 

(4) Define the points E , 0 2 , F and 0 2 as the midpoints of the segments constructed in the previous 
step (3). 

Now the dynamic model is ready and we can play 
with it. Really our model is dynamic (’’alive”) in 
the sense we can dynamically change its parame- 
ters. All initial points D X ,Q and B 2 are moveable - 
they can be moved with the little hand } (with 
the mouse in the PC-Cabri). Playing with this 
model we really see that the third - depending - 
quadrangle is a square (we can convince in it by 
inspection or by measuring the parameters of a 
figure with the Cabri Geometre instruments and/or by applying again the macro on E,F or 0i,0 2 ). Of 
course our experimental assurance needs the deductive proof, it is only the hypothesis, but this hy- 
pothesis is of a great confidentiality. We put aside now the attempts to prove the hypothesis and pro- 
ceed our computer experiments. ( File exp1.92a) 



Experiment 2 (disjointing the vertices) 

Is it really important that the initial squares have a vertex in common? 

Unfortunately, the previous model cannot be modified for these purposes. Therefore we must repeat 
the algorithm discribed above with a little change in the second step - the initial points will be four in- 
dependent points B x , D x , B 2 and D 2 , the diagonal points of the future quadrates. 
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By playing with a new Model 2 we can convince 
ourselves that in this case the resulting quadrangle 
is a square as well. As a result of the experiments 
described we can formulate the next hypothesis: 


Generalization 1 

Two squares A x BiC\Di and A2B2C2D2 with centers O x and 0 2 of identical orientation are given. Let E, 
F, G and //the midpoints of the segments A X A 2 , B X B 2 , CiC 2 and D X D 2 . Then the quadrangle EFGH is a 
square as well. 


I tried to generalize model 1 on the TI-92 and was lucky enough to be successful. You only have to 
define two points B 1 and D 2 on the same spot and then produce the initial squares. Whenever you are 
pointing to that "double point" you are asked "Which object?". You have only to choose one of the two 
answers proposed. The four TI-92 screen shots are showing first only the two squares - first con- 
nected and then disconnected - and in the second row you can see model 1 followed by model 2. 





(files exp2.92a and exp22.92a on the diskette) 


Experiment 3 (moving the middle points) 

Why must the vertices of the resulting figure be the midpoints? 


Maybe the result remains valid in the case of arbitrary points which divide the segments joining corre- 
spondent points in a given proportion? 

Modify our previous model. For that purpose a new macroconstruction divide would be needed, which 
divides a segment in a ratio given by the segment and a point on it. Such a macro can be defined with 
the ideas of the Fale's theorem for example. Then we produce Model 3 using this macro by picking an 
arbitrary point E on the segment A X A 2 and then define all the other points F f G and H applying divide 
as points which divide the correspondent segments in the same ratio. 
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Now again with the mouse (or with the } ) we can change the sizes of the squares and their mutual 
positions and the position of the division points as well. As in the previous case we can easily con- 
vince ourselves that the resulting figure is a square again. Thus we obtain the next generalization. 


Generalization 2 

Two squares A X B X C X D X and A 2 B 2 C 2 D 2 with centers O x and 0 2 of identical orientation are given. Let E, 
F, G and //points of the segments A X A 2 , B X B 2 , CiC 2 and D X D 2 respectively which divide them in the 
same ratio. Then the quadrangle EFGH is a square as well. 



You can find the TI-92 realization on the Diskette of the 
Year 1998 saved as exp3.92a. The macro divide can be 
applied pointing first to the given segment (eg A X A 2 ), then 
to the given dividing point (E) and then to the segment to 
be divided in the given ratio A X E : EA 2 . Josef 
(File exp3.92a) 



Proof in the DERIVE environment 

#1: "SQUARES. MTH, Rakov & Gorokh" 

#2: CaseMode := Sensitive 

Let O x (a,b) and 0 2 (c,d) the centers of the squares A X B X C X D X and ^4 2 5 2 C 2 Z) 2 . n x and n 2 are the vectors 
O x A x and 0 2 A 2 . 

#3: [01 := [a, b] , 02 := [c, d] , nl := [p, q] , ml := [-q, p] , n2 := [r, s] , 

m2 := [ — s, r] ] 

#4: [A1 := 01 + nl, B1 := 01 + ml. Cl := 01 - nl, D1 := 01 - ml] User 

#5: [A2 := 02 + n2, B2 : = 02 + m2 , C2 := 02 - n2, D2 : = 02 - m2 ] User 

Define the function divide: 

X + t • y 

#6: divide (x, y, t) := User 

1 + t 

Declare the coordinates of the points E, F, G and //which divide the segments A X A 2 , B X B 2 , CiC 2 and 
D x D 2 in the ration t : 1 . 

#7: [E := divide (Al, A2, t) , F := divide (Bl, B2, t) , G := divide (Cl, C2, t) , 

H := divide (Dl, D2, t) ] 
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Find the difference of the vectors HE and GF. Simplifying this expression we obtain [0,0]. Conse- 
quently the quadrangle EFGH is a parallelogram. 

#8: (E - H) - (F - G) = [0, 0] User=Simp (User ) 

Evaluate now the scalar product of the vectors GF and FE, which simplifies to 0 . Thus we have 
proved that the quadrangle EFGH is a rectangle. 

#9: (F - G) • (E - F) =0 User=Simp (User ) 

Compare the lengths of the sides GF and FE of this rectangle. For this purpose declare a new function 
SModV(Square of the Module of the Vector). Simplification of expression #1 1 gives 0, consequently 
the rectangle is a square. Note that the square EFGH can degenerate in a point. 

#10: SModV(x) := x • x User 

#11: SModV (F - G) - SModV (E - F) = 0 User=Simp (User ) 


(If you are able to realize Cabri and the CAS of DERIVE on a TI-92 you will proof the conjecture on 
the TI-92. I add the script proof. Loading it in the TextEditor you can execute all lines pressing <t . 
Please note that I had to use the name d i v i d for the function, because divide is in use (for the 
macro) and that I am not allowed to use the reserved system variables c 1 and c 2 . Josef) 
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dot p( x, x) »smodv( x) 
s mo d v ( f ■ g ) - s mo d v ( e ■ f ) 
del var o 1 , o 2 , n 1 , n 2 , ml , m2 
delvar al,a2,bl,b2,ccl,cc2,dl,d2 
delvar divid,smodv 




FEt Y Y FH 

Connand Uiew Execute I 


:rctt 


|s t + c t + a + q -ir -c- crt-D + pn 

L t + 1 t + 1 J 

■e-h-(f-g) [G 0] 

dotP(f - g , e - f ) O 

dotP(x ? x) + snodu(x) Done 

snodu(f - g) - smodu(e - f ) 0 


C : dotp t x , x } -s-snodo t x > 

C : smodu ( f-g > -snodu ( e-f > 
Csdeluar ol ? o2 ? nl ? n2 ? nl ? n2 


Further generalization 

Experiments with arbitrary similar quadrangles allow one to formulate the following generalization: 
Generalization 3 

Let F\ and F 2 two similar figures in the plain of the same orientation (it means that there exists a simi- 
larity / of the first class, which maps the figure F\ onto the figure F 2 . For each point X\ of figure F u 
and a point X 2 =j{X\) of figure F 2 define a point X , which divides the segment X x X 2 in a ratio t with 
t = X\X : XX 2 . Then figure F, which consists of all such points X is similar to the two original figures 
or degenerates to a single point. 


Proof this conjecture! 
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The radius-vector of a point M will be denoted as M in the following. 

I .et A’, and Y\ two arbitrary points of figure F x . LetX 2 - /(AY) and 7 2 - /( 7 ). X and Y are the points di- 
viding the segments XXi and 7 Y 2 in a ratio t. 


Then: 


X = 


X, +tX 2 

i+7 


7 = 


Y l +tY 2 
1 + t 


As a consequence we obtain: 


— XjY[ + iX 2 7 2 
1 + t 


Therefore we have: 



X x Y x 2 +2 tX x Y x -X 2 Y 2 +t 2 X 2 Y 2 

o + n 2 


Denote the angle beween the ray X{Y\ and its image under the similarity / as <p. This angle <p is con- 
stant because the mapping / is a similarity of the first class. With respect to the relation 
| X 2 Y 2 | = k\X l Y l | , where k is the coefficient of the similarity we finally obtain: 


XY 


1 + 2^coscp + t 2 k 2 

(1 + 0 2 


XJi 


2 


Consequently \XY\ = C \X\Yi\, where C is a konstant. 


Thus figure F is similar to figure F\ if k > 0. In the case k= 0 figure F is a point. 

I realized Generalization 3 on the TI-92 working with two 
similar rectangles A 1 B 1 C 1 D 1 and A 2 B 2 C 2 D 2 and the divid- 
ing point E. As you can see the ratio of the lengths of the 
sides is the same for the three figures (= 2.35). Josef 

(File exp4.92a) 
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Remark that in the most general case (Generalization 3) the proof was done without support of the 
computer and was more simple and natural. It is a rather general situation because from the general 
point of view the unimportant details are disappearing and the matter of fact becomes more clear and 
obvious. Nevertheless the computer experiments have played the substantial role in the process of 
generalization. By the way, the proof discussed above could be done on a computer as well but all the 
analytic calculations were so simple that computer support was unnecessary. 


Last remarks 

As it can be clearly seen from the proof given above the fact stays valid in the case of an arbitrary di- 
mension (not only in the case of 2D, but 3D, 4D etc. as well). 

The successfully proved theorem can be used in computer animation for modeling continuous trans- 
formations of figures in plane or space. 

The authors of the article welcome all the enthusiasts of mathematics and computing to continue with 
computer explorations in the discussed problem field in particular in searching for further generaliza- 
tions and applications of the obtained results. 


Literature: [1] Boltyansky V. G. , About one parquet, Mathematics in School, 1984, Nr 1, p65 - 66 
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The squares and TI-NspireCAS (by Josef Bohm) 



This is generalization 2 performed on the Geometry page of Tl-Nspire: 



The ratio can be changed by the slider. Distance x is calculated according the given formula. It is the 
distance between the vertices of the first square to the division points. 
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There are no scripts available with Tl-Nspire. 
tion which can be used in a similar way. 


But instead of this we have the Notes Applica- 


W 


m2: = 


al 

bl 

cl 

dl 

a2 

b2 

c2 

d2 


s 
r 

=ol+nl 

=ol+ml 

=ol-nl 

=o lml 

=o2+n2 

=o2+m2 

=o2-n2 

=o2-m2 


divide (xx^, tt ) : = xx+tt W 
1 +tt 

e:=divide(al,a2,,tt) 

f:=divide(bljb2,tt) 

g:=divide(cljC2,tt) 

h:=divide(dljd2,tt) 

e-h-(f-g) 

dotp(f-ge-f) 

smo dv(xx) : =dotP (xx, xx) 

smo dv(f- g) - smo dv( e -f ) 




e:=divide(al J a2 J ^) 


f:=divide(bl,b2,tt) 


1 +tt 

it l'+a+C' tt+p 
tt + 1 

U- s+Z>+c^ tt+q 
tt+1 


g:=divide(cljC2,tt) 


h:=divide(dljd2,tt) 


-(«■ 

S-Q-C' 

t 

■+-j 


tt+i 


«■ 

r+b+d- 

tt+P 


tt+1 



■o 

1 

ts 

1 

h-L. 

tt+p) 


tt+1 


-(«■ 

1 

1 

tt+q) 


tt+1 


«■ 

s+c+c 

l 

t 

■+-J 


tt+1 


-(« 

1 

1 

• tt+p) 

tt+1 


s-h-(f-g) 


A 


dotP(f-g,e-f) ► OA 
smodv(xx) : =dotP(xx,xx) ► Done 

smodv(f-g)-smodv(e-f) ► OA 


What is to be executed on TI-92 and V200 is here written in so called Math Boxes. 

In the left screen shot all evaluations are suppressed. Now you have the choice evaluating 
one box after the other and following the process step by step 


or 

you can evaluate the “script” from the first row to the last one as a whole. 

I didn’t use variable t because it is used as slider variable in the Geometry App. Same for 
variable x. The “traffic sign” is a warning: Domain of the result might be larger than the do- 
main of the input. 
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3rd Solution: 

Now we will consider the slope of the cross section, i.e. the direction of the tangent and then we will 
calculate an appropriate correction k(x,a ) using the perpendicular distance a: 
g(x)=f{x)+k(x,a), with k(x,a)=a*SQRT(f'(x) A 2 +1). 

Experimenting with various distances a (positive and negative) DERIVE presents terrible surprises 
(outside of the origin parabola). 

F(x) := x a 2/2 


VECTOR ^F 00 + a-jjl + FGoj j, a, , , j VECTOR^FftO + a-J 1 + Ffcojj, a, -6, 6, 0. sj 

4th Solution: 

Let's consider now the fact that we want to plot simultaneously to each point P of the original curve 
the point Q of the parallel curve —> points P{t) and Q(t). Following this idea we find the perfect solu- 
tion of our problem, (only in the Advanced Course). Here we include calculation with vectors. And 
this is my urgent demand: 

Don’t teach Calculus and Linear Algebra separated from each other in upper secondary level but use a 
connected approach where- and whenever it is possible. 

Let v p{t) = [. xp{i ), yp(t)\ the position vector of P(t). Then 
vq(a,t ) = v p{t) + a-vne(t) is the position vector of Q(t) where 
vne(t) is the perpendicular unity vector. 

In our special case we have: 

r\ 2 i ^ 

t a 

VECTOR t, + [-t, 1], a, -B, 3, 0.5 

L 2 J 2 

^ lit + 1 ) j 

DERIVE can show that there are really parallel curves. 

Really for each real al No, there are surprises again - and now inside of the given parabola. But in this 
case we find an obvious reason for that misbehaviour: inside of a narrow parabola there is no space for 
a second one in a distance a. A new point of discussion is born: 
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What is the maximum thickness of the glass? (from now we will only consider parallel curves inside 
of convex curves). 


In generalized form we have the following definitions: 
F(t) := 

VP(t) := [t, F (t) ] 
d 

VR(t) := — VP(t) 
dt 

VN(t) := J- CVR(t))^, (VRCO^j 


AVP(a, t) := VN(t) 

I VN(t) | 

VQCa, t) := VP(t) + AVP(a, t) 

Using a piecewise defined curve - lines and a quarter circle - (instead of the parabola) we can observe, 
that the parallel curve with maximum distance is marking time, even when curve plotting goes on. 

(a = - 1) Therefore we use the ’’plot velocity” of the parallel curve, which is the derivative of vq(a,t). 
If vq\a,t) is the zero vector we can determine the maximum distance b(t ) for each t. 


F (t) 


If t < - J2/2 
t + J2 

If t < 7C1/2) 

7ci - t A 2) 
72 - t 


We return to the parabola from above: 
2 

t 

F(t) := 

2 




By the way, the points on the blue curve (evolute) are the centres of curvature for all parallel curves. 
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F(t) := 

2 

3-t 

+ 1 

2 

Among the DNL#31 -files you will find a Demo-file including all the approaches shown above using 
some other functions. It should inspire for own investigations. Let me close with a question: 

I know: 

The parallel curve of a line is a line, 

The parallel curve of a circle is a (concentric) circle. 

I would like to know: 

Which type of curves allows an equivalent rule? 



Helmut Wunderling 
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You can experiment with other types of curves. 

See below a trig function followed by an exponential function - which was the starting popint of the 
discussion. 



The next two pages show the 2013 approach plotting the parallel curves in the Graph Appli- 
cation of TI-NspireCAS. The curves are produced as loci of a point and then as enveloppe of 
a line parallel to the tangent of the given function graph. 
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The parabola together with the loci of the moveable point. 



Next page gives two screen shots of the exponential “mother curve” and its parallel curves 
produced as loci of a point and as envelope of the parallel to the tangent. 
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The trig function together with its parallel curves produced with TI-NspireCAS. 



The first contributions of the User Forum starting next page are a 1998 discussion about 
SORT-routines. Now these discussions are obsolete because an excellent working SORT 
procedure is built in. Nevertheless the discussion might still be of interest. Josef 
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Reinhard Schaeck, Berlin, Germam 


schaeck@ibm.net 


I have a serious (though trivial) problem with Derive 4.09 for Windows. I want to sort a vector. I hope 
you can reproduce the mysterious result with the following steps: 

1) We load the utility-file NUMBER.MTH (which provides the SORT-function). Test it with 

SORT ( [3,2,5,1,71 ) = [1,2,3,5,71 


which is o.k. 


2) We enter: 

a:= [374, 358, 341, 355, 342, 334, 353, 346, 355,344,349, 330,352,328,336,359, 
361,345,324,386, 335,371,358,328,353,352,366, 354,378,324] 

3) We test the dimension: 

DIMENSION (a) = 30 
which is o.k. 


4) Now we enter SORT(a) 

5) and simplify: 

[648, 656,330,334,335,336, 341,342,344,345,346,349,704,706,354,710, 
716,359,361,366,371,374,378,386] 

which is NOT the sorted vector and has only 24 elements. 

Any idea? Some other vectors with 30 elements work. Sorting should be so simple... 


Jim FitzSimons cherry@neta.com 

That sort utility does not sort duplicate values. 358, 355, 353, 352, 324, 328 are duplicates. The sort 
adds the duplicates 2*324=648, 2*328=656, 2*352=704, 2*353=706, 2*355=710, 2*358=716. 

DERIVE needs a better sort. 


Johann Wiesenbauer, Vienna, Austria J.Wiesenbauer@tuwien.ac.at 

Hi Reinhard, 

You are just ... wait a minute ... the 1274th DERIVE-user who was trapped into thinking that the sort- 
routine in NUMBER.MTH is an all-purpose sort-routine. I assure you (straight from the horse's 
mouth, as it were, since I wrote it) that it isn't one and it was never intended to be one. 

The truth is that it works for distinct nonzero rational numbers only, and in fact, this will cover more 
than 80% of all cases in practice, in particular when it comes to number theory. (Yes, it is true, they 
should have pointed out that restriction in the online help!) 

Rumours are that there will be a built-in SORT-routine in the next release of DERIVE. In the mean- 
time, this is what I can offer you: 

(Rumours were right at this time. In the meanwhile a built-in SORT-routine is available in the 
next versions.) 

SORTO (v) := LIM (TERMS (v-p VECTOR (x_ A (-k_) , k_, v) ) , x_, 1) 

(Sorts any list of distinct nonzero rational numbers.) 

SORT1 (v) := VECTOR (vSUBk_, k_, LIM (DIF (TERMS (EXPAND (VECTOR (x_ A (-v_) , 

v_, 

v) -VECTOR (y_ A k_, k_, 1, DIMENSION (v) ))) , y_) , [x_, y_] , [1, 1])) 

(Sorts any list of nonzero rational numbers.) 
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Sebastiano Cappuccio, Forli 

Connected Mode with the TI-92 


scappucc@spfo.unibo.it 



With DERIVE we can plot a n x 2 matrix to get the 
graphical representation of a set of n points. 

A very useful feature is available for that: the Connected 
Mode, which allows to plot a polygon by plotting the 
matrix of its vertices. 


DERIVE ' s brother, the TI-92 calculator, has not such a feature, but with its powerful programming 
language, we can create it easily. 


Of course, if you put a matrix 3x2, e.g. [[1,1], [3,1], [1,4]], as argument of the program PLOTM, you 
don’t get a triangle; in order to have a triangle you have to “close” the polygon by putting a 4 row in 
the matrix m with the coordinates of the the first point again: [1,1], or you can modify the program. 


Now we can define some new matrix and use matrix multiplication to get, for example, a symmetric 
triangle with respect to the origin or to get a 60° clockwise rotated triangle: 



rFiT'HOY T F2 T FH T 

t t— |Zoom|Trace|ReGraph|Math 


■ plot-pi(pi) 

Done 

r-i 0 ' 
"° S [0 - 1 . 

■ plotri( m ■ os) 

Done 


t\ 




MAIN RA 

AUTD FU 

n: 


P^l^.fcSTcWSi-lprS,I0fcl«r" «_ 1 


f +1/2 



^ + 3^2 

* U2 

■ pi ■ not | y = 

2-J3 + 1/2 



43.,.^ 

-J3 . J 

plotm<ans<l>>l 

MAIN RAD AUTD 

FUNC 



Gunter Scheu, Pfinztal, Germany 

Dear Josef, 

The disussion on page 46 (revised DNL#30) has surprised me very much because "apples are com- 
pared with pears”. How is probability defined there? There is no Laplace experiment and the number 
of outcomes is not finite. 

You know: 

(1) All primes - with exception 2 and 3 - show the form 6n-l or 6n+l (proof is elementar). (Numbers 
of the form (6n- 1, 6n+ 1) are called prime twins, e.g. n= 1, 2.) It is not known if there is an infinite 
number of prime twins. 

(2) Asymptotic estimation of the prime numbers' distribution. 

Jan Vermeylen's question was if the relative frequency tends towards 0.5. Johann Wiesenbauer proved 
this. T. Etchells' remark makes no sense because it refers to a limit. The axioms of probability theory 
(Kolmogoroff) say that 

Pi + P 2 + + p n = 1 which leads in our case to p({2,3}) + p(Z 6n _i) + p(Z 6n +i) = 1. 

If p(Z 6n _i) = 0.5 then the number of primes of the form 6^7-1 would exceed the others of the form 6?/+l. 
Such results cannot be proved, because there is no formula to construct primes, (and following Gold- 
bach's Conjecture there will never be one!). 
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Summary: Be careful to change from finite sets to infinite sets. CAS is at most able to handle finite 
sets. CAS can lead us to conjectures but unfortunately cannot prove them. Best regards to you all, 
Gunter 

DNL: I don't think that Terence's intention was to "prove" the conjecture. In Gunter's sense 
Terence's DERIVE code is affirmating a conjecture. See his comment: "Pure maths wins the 
day!" 


Peter Witthinrich, Luebeck, Germany PWitthinrich@t-online.de 

A comment to "6174 is a Special Number" from DNL # 30 

According [1], a collection of witty brain-twisters, a certain D R Kaprekar from Delavi (India) should 
have discovered 6174 thirty years ago. In [1] you can find the constants from 2 through 7 digits as fol- 
lows: 


k2 

— 

81 

k3 

= 

495 

k4 

= 

6174 

k5 

= 

75933 

k6 

= 

642654 

k7 

— 

7509843 


I have used those ’’Kaprekar's Numbers” several times in teaching Information Technology, because it 
performs as a nice not too large programming problem for Pascal. The natural numbers are represented 
as strings or as arrays and the basic calculations are simulated by appropriate algorithms - like in ele- 
mentary school. In DERIVE it runs much quicker, a fine thing. 

[1] Otto Botsch, In der Werkstatt der Himverzwirner, Aulis-Verlag, Koln 1979, Problem Nr. 79 (page 
115) 


Steven Schonefeld, Auburn, IN, USA schonefelds@alpha.tristate.edu 

Hi Josef, 

It was nice seeing you in Gettysburg. I enjoyed your talk on Dimensional Analysis. I learned 
about this subject in college chemistry class. 

I was impressed by Terence Etchells' Gettysburg talk. One of the things he talked about was the 
implementation of termination tests using DERIVE. It is about time for the procedures in my book [1] 
to be re-written with built-in termination tests. Terence's talk inspired me to do so. Te results of my ef- 
fort are on the enclosed disk. 

Let me say a few words about these procedures. A typical iterative method starts with an initial 
guess, x 0 , and creates a sequence of iterates , {x 0 , xi, x 2 , . . . }, which are computed from the initial 
guess by a function, G(x): xi = G(x 0 ), x 2 = G(xi), etc. 

My feeling is that the calculation of new iterates should be stopped by a two-part termination test. 
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Part one stops the iteration process if two successive x's are close. That is, the iteration stops with x n 
if 


ABS(x n - x n _i) < Stol 

- where Stol is a small, positive number. 

Part two stops the iteration process after a given number of iterations have been performed. That is, 
the iteration stops with x n if 

n > njnax 

- where njnax is a large, positive number. 


Using ITERATES(NE W_V (v),v,v 0 ) : 

With this method, v is a vector containing an iterate x n . As far as I could tell, an iteration counter, 
n, must also be included in the vector v. For example, in my implementation of the Newton-Raphson 
method (2 3 NR.MTH), the function G(x) is 


G(x) := x 


Fjx) 

F'(x) 


and v is of the form [x n , 5 n , F(x n ), n] 


where x n = x n _i + 5 n and 8 n = - F(x n _i) / F '(x n _i). The good news is that the termination test is easy to 
implement - since x n - x n _i = 8 n and n are contained in the current vector. The bad news is that when 
the termination test is satisfied, the last vector v is repeated in the resulting matrix m. In order to avoid 
confusion, I had to Take Off the Last row of m. Hence TOL( ) is the name of the function which does 
this task. 


Using ITERATE(NEW_M(m),m,m 0 ): 

With this method, m is a matrix. Each row of matrix m contains an iterate x n , (usually plus more 
information.) For example, in my implementation of repeated substitution (22REP.MTH), a typical 
row of the matrix, m, is of the form 

[x n , n\ with x n = G(x n _i). 

The good news is that we may determine the iteration counter, n , from DIMENSION(m). (I in- 
cluded n in the typical row of m for pedagogical reasons.) The bad news is that the initial matrix, m 0 , 
must contain at least two iterates for the first part of the termination test 

ABS(m n9 i - m n _i 9 i) < S_tol 

to make sense. Also when using ITERATE( , , ), the function, NEW_M(m), may be quite complicated. 


[1] Numerical Analysis via DERIVE, Steven Schonefeld, 1994, MathWare 


Steven offers an update of his DERIVE-files containing self terminating procedures for 
his book to book owners . If you are interested in the files then please contact Steven 
or me. Thanks a lot for that offer : Josef 
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In revised DNL#20, p 54 you could find "Vieta by Chance", a DERIVE Utility to train solving quad- 
ratics with integer solutions mentally. In my opinion some skills should survive even in times of 
DERIVE and the 77-92, Voyage 200 - and TI-Nspire as well. So I wrote a little 77-program and believe 
it or not, the pupils liked to compete their skill in solving as many quadratics as possible. Josef 

v i e t a ( ) 

P r g m 

Local “ b, x 1 , x 2 , kl, k 2 , vl, v 2 , bz, rz, key 
0 » b z : 0 » r z 
Lbl start 
Cl rl 0 

Disp "Problems: " &s t r i n g ( bz ) &" , correct "&stri ng(rz) 
chart rand( 26) +96) »"b 

Lbl test 

rand( 24) - 12»xl: rand( 34) - 17»x2 
I f x 1 * x 2 =0 or xl= ? x2 
Got o test 

§ x 1 - x 2 » k 1 : xl*x2»k2 
If k 1 <0 Then 
" - " » v 1 
El se 

" + " » v 1 
E n d I f 

If k 2 <0 Then 
" - " » v 2 
Else 

" + " » v 2 

Endl f 

If k 1 =1 Then 
■ " » k 1 
Else 

s t r i n g ( a b s ( k 1 ) ) » k 1 
Endl f 

stri ng(abs(k2))»k2 

Disp expr ( ‘ b&" A 2" &vl&kl&“ b&v2&k2&" = 

0" ) 

Lbl chkl 

Input " 1 . Sol uti on: ", I 1 
Input " 2 . Sol uti on: ", I 2 
If I 1+1 2=xl+x2 and I 1*1 2=xl*x2 Then 
Disp "right" 
r z + 1 » r z 
Else 

Disp "sorry, false "&“b&"l = " &st r i ng( xl) &" , "&“b&"2 = "&stri ng(x2) 

Endl f 


b z + 1 » b z 

Di s p "End = ESC. next = any 
Loop 

get K e y ( ) » k e y 
If k e y =2 6 4 Then 
Goto ende 
Elself k e y V8 Then 
Goto start 
Endl f 
EndLoop 
Lbl ende 
EndPr gm 


- I U.-'l 

■T-- :| . -:,1 

Problems: 1, correct O 


y 2 + 2 y - 15 = O 


1. Solution: 
-5 


2. Solution: 
3 


right 

End = ESC. next = any 


MAIN RAD EKACT 

FUNC 1 3 ■ Ri'l 


Problems: 2, correct 1 

l 2 - 21 ■ 1 + 104 = 0 
1. Solution: 



■ C'k 1 

Problems: O, correct O 


u 2 - b u - 7 = O 


1. Solution: 


MAIN RAD EKACT 

FUNC 


-:|/l 

• ::-k 1 

Problems: O, correct O 


u 2 - b u - 7 = O 


1. Solution: 


1 


2. Solution: 


-7 


MAIN RAD EKACT 

FUNC 


Problems: O, correct O 

u 2 - bu - 7 = O 

1. Solution: 

1 

2. Solution: 

-7 

sorry, false ul=7, u2= - l 
End = ESC. next = any 

MAIN RAD EKACT FUNC 1 E /30 Ijlh'l 


FUNC 1 E /30 
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It is not too difficult to transfer the program to TI-NspireCAS considering the changed form of input 
and output: 

Define vieta()= 

Prgm 

:Local 0b,xxl,xx2,kkl,kk2,vvl,vv2,ss,bz,rz,key,equ 

:bz:=0:rz:=0 

:Lbl start 

:ss:="" 

:0 b:=char(randlnt(97,122)) 

:Lbl vals 

: xx 1 : =randlnt( 1,24)-I2:xx2: =randlnt( 1 , 34)- 1 7 
:If xxl*xx2=0 or xxl+xx2=0: Goto vals 
:kkl:=- xxl-xx2:kk2:=xxl*xx2 
:wl:=when(kkl<0," - + ") 

:vv2:=when(kk2<0," - + ") 

:kkl:=when(kkl=l,"",abs(kkl)) 

:kk2:=abs(kk2) 

:equ:= eb&" A 2"&wl&string(kkl)&0 b&w2&string(kk2)&" = 0" 

:Disp expr(equ) 

:RequestStr "Solutions separated by ,:",ss 
:ss:=expr(" { "&ss&" } ") 

:If ss[l]+ss[2]=xxl+xx2 and ss[l]*ss[2]=xxl*xx2 Then 
:Disp "right" 

:rz:=rz+l 

:Else 

:Disp "Sorry, wrong: "&eb&"l = "&string(xxl)&", "&0b&"2 = "&string(xx2) 

:EndIf 

:bz:=bz+l 

:Disp "Problems: ",bz,", correct: ",rz 
:RequestStr "end = e, next = n: ",key 
:If ord(key)=110:Goto start 
:EndPrgm 


Trrr. i 


RandSeed 11156 


Solutions separated by , 


Solutions separated by 6,-2 
right 

Problems: 1 , correct: 1 
end = e, next = n: n 

e 2 +20* e+99=o| 


m 


Done 





I ° K 1 

Cancel 


I 

H 
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The left column shows a VIETA-session performed 
with the software. 

The right column is a collection of handheld screen 
shots. 





vieta () 


s 2 -4- 5-12=0 

Solutions separated by 6,-2 
right 

Problems: 1 , correct: 1 
end = e, next = n: n 

e 2 +20- e+99=0 

Solutions separated by ,: -11,-9 
right 

Problems: 2 , correct: 2 
end = e, next = n: n 

/ 2 -12-/+32=0 

Solutions separated by ,: 8,4 
right 

Problems: 3 , correct: 3 
end = e, next = n: n 

fr 2 -17- &+30=0 

Solutions separated by ,: -15,-2 
"Sorry, wrong: kl = 2, k2 = 15 " 
Problems: 4 , correct: 3 
end = e, next = n: e 



Another run: 
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Karl-Heinz Keunecke, Kiel, Germany 
kh@keukiel.netzservice.de 

Acceleration due to Gravity ~8.14 m sec ' ? 


1 Introduction 

Data collection by means of computers 
is well known and even in schools a used 
technology. The CBL^ and the CBR^ are 
used in connection with graphic or symbolic 
calculators. These handheld systems are eas- 
ier to handle than PC-based systems. But the 
most important advantage lies in the fact that 
the computer algebra of the 77-92 or the 77- 
89 can be applied to the data. That gives the 
teacher the chance to process the data in 
such a way that matches the curriculum, the 
capability or the interests of his students, his 
own ideas, etc. 

As an example in this paper the gravity 
acceleration is determined from the observa- 
tion of a jumping ball with an ultrasonic 
range meter. The unusual result is explained 
as an effect of the buoancy of the ball in the 
air. 

2 Observation of a jumping ball 
with the CBR 

For our experiments we use a big gym- 
nastics ball with a radius of 34 cm. It jumps 
regularly with a small damping rate and it re- 
flects the sound signal always clearly back to 
the CBR. The ultrasonic ranger was mounted 
at a height of 2.2 m. The ball is held at a dis- 
tance of 50 cm from the CBR. As soon as the 
measurement starts it is released and falls 
down to the ground. 


The distance from ground is shown in 



Figure 1 : Height of a jumping ball 

5 rebounds at the ground can be recognized. 

In order to inspect the physical process 
more detailed, a section of figure 1 is con- 



Figure 2: Balls height as a function of time 

The calculated velocity and acceleration are 
displayed in figures 3 and 4. 



Figure 3: Velocity of the jumping ball 








Figure 4: Acceleration of the jumping ball 


Figure 5: Quadratic regression function 


From figure 3 it can be seen that the ab- 
solute value of the velocity increases linearly 
to the ground. There the ball is stopped 
within a short period and is accelerated then 
in opposite direction. The ball moves up 
again with linear velocity until it arrives at 
maximum height where the velocity is zero. 
From there the described ball drop begins a 
second time. The acceleration in figure 4 in- 
creases to about 20g during the ball is re- 
bounded from the ground. It can also be seen 
that the rebounce produces oscillations of the 
ball for about half a second. 



Figure 6: Distance-time diagram (squares) 
and quadratic regression function 

Using the same procedure to find a lin- 
ear regression for the velocity the equation 

v = 8.14 t . (2) 


3 Data processing with the TI-92 

After the observations have been dis- 
cussed by the students, they dealt with the 
question: What functions describe the height 
and the velocity of the ball? 

These questions can easily be solved by 
using the tools of statistical calculations of 
the TI-92. For that purpose the lists 11 to 14 
are read into a table by means of the 
Data/Matrix Editor. Then a regression func- 
tion can be calculated. Choosing a quadratic 
regression the calculated function is shown 
in figure 5. The coefficients b and c of the 
quadratic term are negligible and therefore 
the describing function is approximately: 


is yielded. The data and the regression func- 
tion are plotted in figures 6 and 7. 



Figure 7: Velocity-time diagram (squares) 
and the regression line 

The regression function for the distance and 
for the velocity as well yields an acceleration 
of 8.14 m sec' 2 . 
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Now the students are faced to the ques- 
tion: Why is the gravity acceleration clearly 

_2 

smaller than 9.81 m sec ' ? 

Friction cannot be the reason. The influ- 
ence of friction increases when the velocity 
goes up. Therefore it would reduce the ac- 
celeration at higher speed. This cannot be 
seen from the data for the velocity increases 
linearly through the period of observation. 

But the buyoant force of the ball in the air 
reduces the gravity force to a certain amount 
which can be determined by using the prin- 
ciple of Archimedes. 

The reduced gravity force of the ball in 
air G' equals the difference of the gravity 
force in vacuum G and the buoyancy of the 
ball A : 

G' = G - A (3) 

The ball is buoyed up with a force A equal to 
the gravity force of the displaced air 

A = m d i S g (4) 


Now the reduced gravity acceleration g' can 
be determined from equation (6): 


g =g 


m + m ■ 

env air 


The mass m air can be calculated from the ra- 
dius r of the ball and the density p mr of the 
air, 

77 471 3 

m air = m dis = Pair V = Pair ~ r , 


r= 34cm, p air = 1.29- 


m air= 0.212kg 


The mass of the ball determined by a pair of 
scales, equals under the made assumptions 
the mass of the envelope, 

m env = 1.158 kg. 

With these results the reduced gravity accel- 
eration can be calculated from equation (8): 


m | 1.158 kg 

sec 7 ll.l58kg + 0.212 kg 


(4) 

g' = 9.8 1 • 

m 

2 

= mg’ 


sec 

i is the 

= 8.29 

m 


of the 


sec 2 


ball and the mass m a i r of the air included in 
the ball: 

m Hlenv "F Hlair (5) 

Then equation (3) yields: 

( m c nv "F m a ii -) g g “1“ m a ir g Wilis g (6) 

The pressure in the ball is not much above 
the atmospheric pressure. Therefore it can be 
assumed that the mass of the air in the ball 
equals nearly the mass of the air which is 
displaced by the ball, 

[1] Calculator Based Laboratory, Texas Instruments 

[2] Calculator Based Ranger, Texas Instruments 


The calculated acceleration does not 
agree totally with the observed acceleration. 
The reason probably is the assumption of 
equation (4). The mass of the air included in 
the ball is, because of the pressure in the 
ball, greater than the air being displaced by 
the ball. 






